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Ionospheric  Effects  on  Earth -Space  Propagation 


1.  INTRODUCTION 

Radio  waves  undergo  several  effects  when  they  pass  through  the  earth's  ionos¬ 
phere.  One  of  the  most  important  effects  that  the  ionosphere  has  upon  radio  waves 
that  traverse  it  is  a  retardation,  or  group  delay,  on  the  modulation  or  information 
carried  on  the  radio  wave  that  is  due  to  its  encounter  with  the  free,  thermal 
electrons  in  the  earth’s  ionosphere.  Other  effects  the  ionosphere  has  on  radio  waves 
include:  (1)  RF  carrier  phase  advance;  (2)  Doppler  shift  of  the  RF  carrier  of  the 
radio  wave;  (3)  Faraday  rotation  of  the  plane  of  polarization  of  linearly  polarized 
waves;  (4)  angular  refraction  or  bending  of  the  radio  wave  path  as  it  travels  through 
the  ionosphere;  (5)  distortion  of  the  waveform  of  transmitted  pulses,  and  (6)  ampli¬ 
tude  and  phase  scintillation.  With  the  exception  of  scintillation  effects,  which  are 
not  described  in  this  report,  all  the  other  effects  listed  here  are  proportional,  at 
least  to  first  order,  to  the  total  number  of  electrons  encountered  by  the  wave  on  its 
passage  through  the  ionosphere  or  to  their  time  rate  of  change.  In  fact,  phase  scin¬ 
tillation  also  is  merely  the  short  term  time  rate  of  change  of  total  electron  content 
(TEC)  after  the  longer  term  variations  have  been  removed. 

In  this  report,  a  short  description  is  given  of  each  ionospheric  TEC  effect  upon 
radio  waves,  along  with  a  representative  value  of  the  magnitude  of  each  of  these 
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effects  under  normal  ionospheric  conditions.  This  is  followed  by  a  discussion  of  the 
important  characteristics  of  average  ionospheric  TEC  behavior  and  the  temporal 
and  spatial  variability  of  TEC. 

1.1  Group  Path  Delay 

The  additional  time  delay,  over  the  free  space  transit  time,  of  a  signal  trans¬ 
mitted  from  above  the  ionosphere  to  a  user  on,  or  near,  the  earth’s  surface  is  given 
by  Millman:  ^ 

4  0  s 

At  =  — ‘nr-  TEC  (seconds)  tl) 

cf 

where  the  TEC  is  the  total  number  of  electrons  along  the  path  from  the  transmitter 

to  the  receiver,  c  is  the  velocity  of  light  in  m/sec,  and  f  is  the  system  operating 

frequency  in  Hertz.  The  TEC  is  generally  expressed  as  the  number  of  electrons  in 

a  unit  cross  section  column  of  one  square  meter  area  along  this  path. 

A  plot  of  time  delay  vs  system  operating  frequency  for  TEC  values  from  101  to 
10  *3 

10  el/m  is  given  in  Figure  1.  These  two  values  represent  the  extremes  of  ob¬ 
served  TEC  in  the  earth’s  ionosphere.  Note  that,  at  a  system  operating  frequency 

1 8 

of  1  GHz,  for  example,  a  TEC  of  10  ,  a  value  frequently  exceeded  in  many  parts 

of  the  world,  would  produce  a  time  delay  of  134  nanoseconds  or  132  ft  (40.2  m)  of 
range  error.  At  a  system  operating  frequency  of  100  MHz  this  same  TEC  value 
would  produce  a  range  error  of  over  13,  000  ft  or  4  km!  Obviously,  the  TEC  param¬ 
eter  is  of  potentially  great  importance  to  precision  satellite  ranging  systems. 

1.1.1  TWO-FREQUENCY  IONOSPHERIC  TIME 
DELAY  CORRECTIONS 

If  the  navigation  or  ranging  system  bandwidth  is  large  enough  so  that  two,  fairly 
widely  spaced  bands  can  be  used  for  ranging,  the  ionospheric  time  delay  error  can 
be  reduced  to  an  acceptable  level  automatically,  and  can  be  made  tronsparant  to  the 
system  user.  Because  the  ionospheric  delay  is  a  function  of  frequency  we  can  write; 

At.  =  TEC  ,  A t„  =  -V  ”  TEC 

cf? 

where  Atj  is  the  Ionospheric  error  on  frequency  f^  and  At^  is  the  ionospheric  error 
on  the  frequency  fg.  If  the  normal  system  operational  frequency  is  and  we  chtn>se 
f^  at  a  lower  frequency  for  ionospheric  correction  purposes,  we  obtain 


t 


See  Reference  7, 


page  31'.  1 


(1 


Figure  1.  Time  Delay  vs  Frequency  for  Various 
Values  of  TEC 


0(At)  =  -  '  TEC  <l/f9  -  1/f2,) 
c  z  \ 


or 


Atl  =  f2/(fl  '  f2]  6(At>  1  (2) 

The  value  6(At)  is  obtained  from  the  difference  of  the  simultaneous  measurements 
of  the  total  range,  including  ionospheric  time  delay,  at  the  two  frequencies,  f  and 
f9,  since  the  geometric  distance  is,  of  course,  the  same  at  all  frequencies.  The 
quantity  f 0 / (f ^  -  f9)  is  called  the  ionospheric  scaling  factor.  For  small  ratios  of 
fj/f^  this  factor  is  much  larger  than  unity  and  the  required  precision  of  the  differ¬ 
ential  measurement  may  be  unreasonably  large.  A  plot  of  this  quantity,  normalized 
by  fj,  is  given  in  Figure  2.  In  this  derivation  the  contribution  of  receiver  noise  to 
the  differential  measurement  accuracy  has  not  been  considered. 
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Figure  2.  Ionospheric  Seating  Factor  vs  Ratio  of 
Secondary  (lower)  to  Primary  (higher)  Frequency 


1.  1.2  AN  EXAMPLE  OF  A  TWO-F  REQUENCY 
IONOSPHERIC  TIME  DELAY  SYSTEM 

The  Department  of  Defense  is  currently  testing  an  advanced  navigation  system, 

12 

called  the  NAVSTAR-Global  Positioning  System  (GPS),  ’  which  uses  coherently 
derived,  identical  modulation  on  two  carrier  frequencies,  called  LI  and  L2,  to 
measure  the  ionospheric  group  path  delay  directly  and  thereby  correct  for  ionos¬ 
pheric  time  delay.  The  ratio  of  frequencies  used  in  the  GPS  system  is  exactly 
154/120,  with  the  higher  frequency  (LI)  at  1575  MHz.  The  two  carrier  frequencies 
transmitted  by  the  GPS  system  are  the  154th  and  120th  harmonics  of  10.  23  MHz. 

This  10  MHz  frequency  is  bi-phase  modulated  on  both  carriers  with  a  psuedo  random 

2 

code  resulting  in  a  [(sin  x)/x]  shaped  spectrum  of  width  20  MHz  to  the  first  nulls. 

A  user  with  knowledge  of  the  transmitted  code  collapses  the  received  spectrum  to 
equivalent  carriers  with  10  MHz  modulation.  The  10  MHz  modulation  is  transmitted 
with  a  known  phase  difference  on  the  two  carriers,  and  the  received  modulation  phase 
difference  is  a  direct  measure  of  the  ionospheric  group  path  delay. 


1.  Demaro,  R.  P.  (1981)  The  all-purpose  satellite,  IEEE  Spectrum.  Hh35-40, 

2.  Milliken,  R.  J. ,  and  Zaller.  C.  J.  ( 1978)  Principle  of  operation  of  NAYSTA  R 

and  system  characteristics,  Navigation:  Journal  of  the  Institute  of  Navigation, 
25 (No”.  2 ) : 95 -  1 06 . 
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l  or  the  GPS  carrier  and  modulation  frequencies  the  ionospheric  group  path 
delay  at  frequency  LI,  as  obtained  from  Eq.  (2)  is: 

Att  =  -1.  5457  <5(At) 

where  5(At)  is  the  difference  between  the  ionospheric  ti-me  delay  measured  at  the 
two  frequencies.  This  difference  in  range  is  directly  related  to  absolute  ionospheric 
time  delay  as,  of  course,  the  satellite  is  at  the  same  range  at  both  frequencies.  The 
only  frequency  dependent  parameter  in  range  measurements  is  the  ionospheric  time 
delay  effect,  assuming  the  satellite  transmitted  modulation  phase  at  LI  and  L.2  is 
known  and  the  receiving  system  frequency  dispersive  characteristics  can  be  inde¬ 
pendently  measured  and  corrected  for. 

18  2 

For  a  typical  daytime  high  solar  activity  TEC  value  of  10  el/m  column  the 
5(At)  measured  by  a  GPS  receiver  would  be  54  nanoseconds  or  53.  1  ft  (16.  2  m)  of 
ionospheric  error.  For  a  direct  measure  of  absolute  TEC  from  the  modulation 
phase  delay  at  L2  minus  L.l  we  have: 

TEC  =  1.  845  *  10 16  ■:=  6(At) 

where  6(At)  is  measured  in  nanoseconds. 

Since,  at  10.  23  MHz,  one  complete  cycle  of  modulation  phase  of  360°  is  97.  75 

16  2 

nanoseconds,  we  obtain:  TEC  =  0.7745  X  10  el/m  per  degree  of  10.  23  MHz 

16  2 

modulation  phase  difference,  or  TEC  =  278.8  X  10  el/m  per  cycle  of  modulation 
phase  difference.  Thus,  the  cycle  ambiguity  in  absolute  values  of  TEC  is  trivial  to 
resolve  using  the  GPS  system  as  a  means  of  determining  ionospheric  time  delay. 

Absolute  ionospheric  time  delay  measurements  can  be  made  with  an  accuracy 
better  than  1  to  2  nanoseconds,  depending  upon  the  received  signal-to-noise  ratio 
on  both  frequencies.  For  the  power  levels  transmitted  by  the  GPS  satellites,  an 
omni-directional  receiving  antenna,  and  a  receiver  with  a  modulation  tracking  band¬ 
width  of  approximately  15  Hz,  the  differential  modulation  phase  has  been  measured 
to  within  approximately  ±  2  nanoseconds.  The  contribution  of  receiver  noise  for  the 
two-frequency  ionospheric  time  delay  corrections  on  the  GPS  system  has  been  con- 

3 

sidered  by  Cretcher. 

3.  Cretcher,  C.  K.  (1975)  Ionospheric  effects  in  NAVSTAR-GPS,  in  Effects  of  the 
Ionosphere  on  Space  Systems  and  Communications,  J.  M.  Goodman,  Ed.  , 

Naval  Research  Laboratory,  Washington,  D.  C.  US  Government  Printing 
Office,  Stock  No.  008-051-00064-0. 
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1.2  RF  Carrier  Phase  Ad%-ance 

In  addition  to  group  path  delav,  or  modulation  time  delay,  over  the  free  space 
delav,  the  phase  of  the  carrier  of  radio  frequency  transmissions  is  changed  by  the 
ionosphere.  The  RF  phase  is  advanced  with  respect  to  its  phase  in  the  absence  of 
an  ionosphere.  This  effect  is  extremely  important  in  determining  space  object 
velocities  bv  means  of  range  rate  measurements.  The  amount  of  phase  increase 
or  phase  path  decrease  can  be  expressed  as: 


AT 


1.  34  X  10 
- r - 


TEC 


(cvcles) 


(3) 


2 

where  f  is  the  system  operating  frequency  in  Hertz,  and  TEC  is  in  el/m  column. 

In  practice,  the  amount  of  this  phase  advance  cannot  readily  be  measured  on  a 
single  frequency  and  two,  coherently  derived,  frequencies  are  required  for  this 
measurement. 

1.2.  1  DIFFERENTIAL  CARRIER  PHASE 

In  addition  to  the  dual  frequency  identical  modulation  transmitted  from  the  GPS 
satellites  for  ionospheric  group  path  correction,  these  satellites  also  transmit  two. 
coherently -derived  carrier  frequencies  for  ionospheric  differential  carrier  phase 
measurements.  For  the  pair  of  frequencies  used  by  GPS,  approximately  1.2  and 
1.6  GHz.  the  differential  phase  shift,  referenced  to  the  lower  frequency,  is: 

-7  2 

AO  =  *'  34  *  10 _ x  — — — tt— -  X  TEC  (cycles)  (4) 

UT 

where  m  =•  f  / F  =  1.2833.  For  the  GPS  frequencies  Ab  =  4.31  X  10  17  *  TEC  or 
16  ^  ^  2 

2.32*  10  el/m  per  complete  2 ir  cycle  of  differential  carrier  phase  between  LI 
and  L2.  measured  at  L2.  The  differential  carrier  phase  [Eq.  (4)]  is  related  to  the 
differential  modulation  phase,  [Eq.  (2)}  simply  by  the  ratio  of  carrier  to  modulation 
frequencies.  With  a  reasonable  carrier  signal -to-noise  ratio,  this  differential 
carrier  phase  can  be  measured  to  within  a  few  degrees,  or  less  than  approximately 
0.  04*  1016  el/m2.  Since  the  TEC  is  generally  much  greater  than  2.  32  X  1016, 
corresponding  to  2  v  of  differential  carrier  phase,  there  is  a  2  nrr  ambiguity  in  the 
differential  phase  measurement. 

The  differential  carrier  method  of  measuring  TEC  cannot,  in  practice,  be  used 
to  measure  absolute  values  of  TEC  by  itself  due  to  the  large  2  nn  phase  ambiguity 
in  the  measurement,  but  this  is  not  important  for  navigation  systems  which  require 
a  correction  only  for  range  rate  errors  due  to  the  ionosphere  between  two  measure¬ 
ment  times. 
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The  US  Navy  Navigation  Satellite  System,  NNSS.  determines  position  for 
stationary  and  slowly  moving  vehicles  by  measuring  satellite  transmitted  RF 
carrier  phase  changes  as  a  function  of  low -orbit  satellite  motion  across  the  sky. 
This  method  of  positioning  requires  only  range  rate  information.  The  primary 
NNSS  frequency  is  400  MHz.  A  second  RF  carrier  at  150  MHz  is  used  only  for 
ionospheric  range  rate  corrections.  While  various  techniques  have  been  proposed 
for  determining  the  absolute  TEC  from  the  differential  carrier  phase  information 
received  from  the  NNSS  satellites,  they  all  involve  assumptions  concerning  some 
a-priori  knowledge  of  the  ionosphere,  and  they  cannot  be  used  in  the  general  case. 

As  an  ionospheric  monitoring  tool  the  combination  of  differential  carrier  phase 
and  differential  modulation  phase  provides  an  excellent  means  of  determining 
ionospheric  electron  content  along  the  ray  path  to  the  satellite.  The  absolute  value 
of  TEC  can  be  determined  by  the  group  delay  technique  and  relative  TEC  changes 
can  be  measured  with  great  accuracy  by  the  differential  carrier  phase  technique. 

1.  2.  2  SECOND  DIFFERENCE  OF  CARRIER  PHASE 

The  second  difference  in  phase  between  an  RF  carrier  and  that  of  its  upper  and 
lower  sidebands  can  be  used  to  measure  absolute  values  of  TEC,  as  described  by 

4 

Burns  and  Fremouw.  If  three,  coherently  derived  frequencies,  f-fm*  f,  and 

f-*-  f  are  transmitted  the  second  difference  of  phase  is  given  by: 
m 


A„  O  =  (4 


-0  ) 
c 


-  (0. 


-  V 


0+0. 

u  I 


2X0 


From  Eq.  (3), 


AO  =  1 '  34  *-1Q  -  X  TEC  (cycles) 


thus 


a2o  - 


2.  68  X  10"7  f2 

2  TEC  (cycles) 


nr  -  r  > 

m 


when  f2  >  >  f  2 
m 


A20  ■ 


2.  68  X  10‘7  f2 

- —  TEC  (cycles) 


f 


T 


(5) 


*See  References  30  and  31,  Reference  page  29. 

4.  Burns,  A. A.,  and  Fremouw,  E.  J,  ( 1970)  A  real -time  correction  technique 
for  transionospheric  ranging  error,  IEEE  Trans,  on  Antennas  and  Propag, 
AP-18(No.  6): 
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For  a  carrier  frequency  of  100  MHz  a  modulation  frequency  of  1.  93  MHz  would 

18  2 

be  required  to  give  2  77  of  second  differential  phase  for  a  TEC  value  of  10 1  el/m  . 

18  2 

A  value  of  A  2  q  of  2  k  for  10  el/m  is  a  reasonable  compromise  between  the 
requirement  for  minimizing  chances  of  an  ambiguity  in  absolute  TEC  and  accuracy 
in  measuring  TEC  relative  changes.  The  second  difference  of  carrier  phase  has 
been  used  with  the  DNA-002  satellite  to  make  estimates  of  the  absolute  values  of 
TEC.  ° 

1.3  Doppler  Shift 

Since  frequency  is  simply  the  time  derivative  of  phase,  an  additional  contribu¬ 
tion  to  geometric  Doppler  shift  results  due  to  changing  TEC.  This  additional  fre¬ 
quency  shift  is  generally  small  compared  to  the  normal  geometric  Doppler  shift, 
but  can  be  computed  bv: 

Af  =  . y-1---  TEC  (Hertz) 

For  high  orbit  satellites  where  the  diurnal  changes  in  TEC  are  greater  than 
geometric  ones,  an  upper  limit  to  the  rate  of  change  of  TEC  is  approximately 

16  9 

0.  1'  10  el/m“  per  second.  This  value  yields  an  additional  frequency  shift  of  less 
than  one-tenth  of  a  Hertz  at  1.6  GHz  which  would  not  be  significant  compared  with 
a  typical  required  receiver  loop  bandwidth  of  at  least  a  few  Hertz.  At  400  MHz  a 
similar  rate  of  change  of  TEC  would  produce  a  frequency  shift  of  approximately 
0.  3  Hz,  probably  still  not  significant. 

During  times  of  severe  phase  scintillation,  which  can  occur  even  at  GHZ  fre¬ 
quencies,  the  TEC  may  not  change  in  a  consistent,  rapid  manner  to  yield  greater 
ionospheric  Doppler  shifts,  but  the  phase  of  the  incoming  RF  signal  can  have  a  large 
random  fluctuation  superimposed  upon  the  changes  associated  with  the  normal  rate 
of  change  in  TEC.  This  large,  random  component  may  actually  spread  out  the 
spectrum  of  the  received  signal  sufficiently  to  cause  the  receiver  to  lose  phase  lock, 
as  the  receiver  signal  phase  may  have  little  energy  remaining  in  the  carrier,  and 
instead  may  be  spread  over  several  Hertz,  with  little  recognizable  carrier  remain¬ 
ing.  A  knowledge  of  phase  scintillation  rates  is  required  to  determine  the  spread 
of  received  signal  phase. 

5.  Fremouw,  E.  J. ,  Leadabrand,  R.  L. ,  Livingston,  R.  C. ,  Cousins,  M.  D.  , 

Rino,  C.  L.  ,  Fair,  B.  C. ,  and  Long,  R.A.  (1978)  Early  results  from  the 
DNA  wideband  satellite  experiment-complex  signal  scintillation.  Radio  Science, 
13(l\o.  1):  167-187. 
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1.4  Faraday  Polarization  Rotation 


When  a  linearly  polarized  radio  wave  traverses  the  ionosphere  the  wave  under¬ 
goes  rotation  of  the  plane  of  polarization.  At  frequencies  of  approximately  100  MHz 
and  higher  the  amount  of  this  polarization  rotation  can  be  described  by: 


2. 36  X  10 


-5 


/■ 


B  cos  9  N  dl 


(6) 


where  the  quantity  inside  the  integral  is  the  product  of  electron  density  times  the 
longitudinal  component  of  the  earth's  magnetic  field,  integrated  along  the  radio 
wave  path.  Many  ionospheric  workers  have  used  this  effect,  named  for  Michael 
Faraday  who  first  observed  polarization  changes  in  an  optical  experiment,  to  make 
measurements  of  the  TEC  of  the  ionosphere.  Since  the  longitudinal  magnetic  field 
intensity  changes  much  slower  with  height  than  the  electron  density  of  the  ionosphere, 
the  equation  can  be  rewritten  as: 

*  -%r  Bl  *  TEC  (7) 


where  B.  =  B  cos  9  is  taken  at  a  mean  ionospheric  height,  usually  near  400  km, 

K  =  2. 36  X  10  and  TEC  is  /  N  dl.  Typical  values  of  polarization  rotation  for 
northern  mid-latitude  stations  viewing  a  geostationary  satellite  near  their  station 
meridian  are  given  in  Figure  3  as  a  function  of  system  frequency  and  total  electron 
content.  In  fact,  the  largest  portion  of  TEC  data  available  today  from  stations 
throughout  the  world  have  come  from  Faraday  rotation  measurements  from  geo¬ 
stationary  satellite  VHF  signals  of  opportunity. 

For  satellite  navigation  and  communication  designers,  however,  the  Faraday 
polarization  rotation  effect  is  a  nuisance.  If  a  linearly  polarized  wave  is  trans¬ 
mitted  from  a  satellite  to  an  observer,  on,  or  near,  the  surface  of  the  earth,  the 
amount  of  polarization  rotation  may  be  nearly  an  odd  integral  multiple  of  90  degrees, 
thereby  giving  no  signal  on  the  receiver's  linearly  polarized  antenna,  unless  the 
user  is  careful  to  realign  his  antenna  polarization  for  maximum  received  signal. 

As  shown  in  Figure  3,  at  4  GHz,  a  commercial  satellite  transponder  frequency 
band,  the  amount  of  Faraday  rotation  can  be  a  tenth  of  a  radian,  well  in  excess  of 
that  required  for  dual,  linear  orthogonal  channel  separation. 

The  Faraday  rotation  problem  is  overcome  by  the  use  of  circular  polarization 
of  the  correct  sense  at  both  the  satellite  and  at  the  user's  receiver.  Generally  the 
mobile  user  finds  it  difficult  to  utilize  circular  polarization  due  to  the  continual 
vehicle  directional  changes;  thus  he  settles  for  a  received  linear  polarization.  The 
3  dB  loss  between  transmitted  circular  polarization  and  receiver  linear  polarization 
is  a  necessary  price  to  pay  for  user  antenna  maneuverability  and  simplicity. 


Figure  3.  Faraday  Polarization  Rotation  vs 
Frequency  for  Various  Values  of  TEC 


1.5  Angular  Refraction 


The  refractive  index  of  the  earth's  ionosphere  is  responsible  for  the  bending  of 

radio  waves  from  a  straight  line  geometric  path  between  satellite  and  ground.  This 

angular  refraction  or  bending  produces  an  apparent  higher  elevation  angle  than  the 

6 

geometric  elevation.  Millman  and  Reinsmith  have  derived  expressions  relating  the 

refraction  to  the  resultant  angular  bending.  Perhaps  the  easiest  expressions  to 

6 

use,  as  given  by  Millman  and  Reinsmith  relate  the  ionospheric  range  error  to 
angular  refraction: 


AE 


R  +  r  sin  E  ( r  cos  E  ) 
o  o  o  o 

- - - - 

h.  (2  r  +  h.)  +  r  sin  E 
i  o  i  o  o 


AR 

“TT 


(8) 


where  E  is  the  apparent  elevation  angle,  R  is  the  apparent  range,  AR  is  computed 
o  A 

from  A  R  =  (40.  3/f  )  *  TEC,  rQ  is  the  earth's  radius,  and  h^  is  the  height  of  the 
centroid  of  the  TEC  distribution,  generally  between  300  and  400  km. 


G.  Millman,  G.  IT. ,  and  Reinsmith,  G.M.  (1974)  An  Analysis  of  the  Incoherent 
Scatter-Faraday  Rotation  Technique  for  Ionospheric  Propagation  Error 
Correction,  Report  No.  RV4^Ml(2,  General  Electric  Company, 

Syracuse,  New  York. 
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For  low  elevation  angles  and  satellites  well  above  most  of  the  ionization, 

R  >  >  r  sin  Eq,  and  the  angular  refraction  can  be  expressed  as: 

cos  E 

AE  =  -■$-h  AR  .  (9) 

1 
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Typical  values  of  elevation  refraction  error  for  a  TEC  of  10  el/m  column 
are  shown  in  Figure  4  for  several  frequencies.  Note  that,  at  the  lowest  frequency, 
100  MHz,  near  the  horizon  the  refraction  is  well  over  1.5  degrees!  The  curves 
shown  in  Figure  4  have  been  constructed  using  the  approximation  derived  by  Millman 

g 

and  Reinsmith  for  low  elevation  angles  given  in  Eq.  (9). 
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Figure  4.  Refraction  in  Elevation  Angle  vs  Elevation  Angle  for  Indicated 
Frequencies  and  Values  of  TEC 


Generally,  the  range  error  Itself  is  the  main  ionospheric  problem  for  advanced 
navigation  systems,  and  elevation  angle  errors  are  insignificant.  Satellite  detection 
radar  systems,  on  the  other  hand,  do  have  the  requirement  to  know  accurate  pointing 
elevation  angles  for  their  large  aperture  arrnys,  though  generally  the  accurate 
tracking  is  done  by  using  range  rate  information,  and  elevation  angle  is  of  secondary 
importance  as  long  as  the  beamwldth  of  the  antenna  is  large  enough  to  see  the  target. 
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Errors  in  the  azimuth  of  radio  waves  transmitted  through  the  ionosphere  can 
also  occur;  they  depend  upon  azimuthal  gradients  in  TEC  which  are  generally  small 
and  which  usually  can  be  neglected  in  practical  cases. 


1.6  Distortion  of  Pulse  Waveforms 


Two  characteristics  of  the  ionosphere  can  produce  distortion  of  pulses  of  UF 

energy  propagated  through  it.  The  dispersion,  or  differential  time  delay  due  to 

7  2 

the  normal  ionosphere,  as  derived  by  Millman  is  proportional  to  1/f  ,  and  pro¬ 
duces  a  difference  in  pulse  arrival  time  across  a  bandwidth  Af  of: 


At 


80.6  *  10 
3 - 


c  f 


Af  X  TEC 


(10) 


where  c  is  the  velocity  of  light  in  m/sec,  f  and  Af  are  expressed  in  Hertz,  and  TEC 
2 

is  in  el/m  column.  The  dispersive  term  for  pulse  distortion  is  thus  proportional 
to  TEC.  When  the  difference  in  group  delay  time  across  the  bandwidth  of  the  pulse 
is  the  same  magnitude  as  the  width  of  the  pulse  it  will-be  significantly  disturbed  by 

g 

the  ionosphere.  Millman  and  Bell  also  derived  mathematical  relationships  for 
ionospheric  dispersive  effects  on  an  FM  Gaussian  shaped  pulse. 

In  addition  to  pulse  distortion  by  the  dispersive  effects  due  to  the  TEC  of  the 
normal  background  ionosphere,  radio  pulses  are  also  modified  by  scattering  from 

9 

ionospheric  irregularities.  Yeh  and  Liu  have  computed  pulse  mean  arrival  time 
and  mean  pulsewidth  due  to  both  dispersion  and  scattering. 


2.  IONOSPHERIC  TOTAL  ELECTRON  CONTENT  (TEC) 


2.1  Average  TEC  Behavior 

The  ionospheric  parameter  responsible  for  the  effects  described  in  Section  1 
of  this  report  is  the  total  number  of  free  electrons,  TEC,  or  its  rate  of  change, 
along  the  path  from  a  satellite  to  a  ground  station.  The  greatest  contribution  to 
TEC  comes  from  the  F2  region  of  the  ionosphere.  A  typical  daytime  mid-latitude, 
high  solar  maximum  electron  density  profile  is  illustrated  in  Figure  5.  The  curve 
on  the  left  side  of  Figure  5  is  the  log  of  Ng  plotted  vs  height  as  normally  shown  by 

7.  Millman,  G.  H.  (1965)  in  Modern  Radar,  Ed.  by  R.  S.  Berkowitz,  John  Wiley  & 

Sons,  New  York. 

8.  Millman,  G.  H. ,  and  Bell,  C.  D.  (1971)  Ionospheric  dispersion  of  an  FM  electro¬ 

magnetic  pulse,  IEEE  Trans.  Antennas  Propag.  AP  -  19(No.  1). 

9.  Yeh,  K.C. ,  and  Liu,  C.  II.  ( 1979)  Ionospheric  effects  on  radio  communication  and 

ranging  pulses,  IEEE  Trans.  Antennas  Propag.  AP-27(No.  6). 
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«l/cm*  el/cm5  unit*  of  I06 


Figure  5.  Typical  Profile  of  Electron  Density  vs  Height.  In  (a),  log  Ne  is 
plotted;  in  (b)  is  plotted  on  a  linear  scale 


ionospheric  workers.  Since  the  TEC  is  represented  by  fhe  area  under  the  curve  of 

linear  plot  of  Ne  vs  height,  the  right  hand  plot  of  Figure  5  illustrates  the  actual 

linear  plot.  Note  that  most  of  the  contribution  to  TEC  occurs  near  the  peak  of  the 

F2  region,  which  is  fortunate,  as  ground-based  ionosondes  have  been  used  since  the 

1930's  to  make  continuous,  routine  measurements  of  the  density  at  the  peak  of  the 

F2  region.  Ionosondes  measure  foF2,  which  is  related  to  N  by: 

max 

(foF2)2  =  80.6  N 

max 

g 

where  foF2  is  in  MHz,  and  N  is  in  units  of  10  el/cc. 

max 

In  the  1950's  and  1960’s  continuing  to  a  more  limited  extent  even  today,  upwards 
of  150  ionosondes  were  operated  to  provide  improved  prediction  capability  for  long 
distance  high  frequency  propagation  by  means  of  ionospheric  refraction.  Various 
models  of  foF2  were  developed  for  this  purpose,  one  of  the  more  popular  ones  being 


17 


commonly  known  as  ITS-7  8^  after  the  report  number  which  described  the  model. 
This  model,  among  other  things,  characterized  the  10-day  average  world-wide 
behavior  of  foF2  by  Fourier  temporal  components  and  Legendre  polynomial  geo¬ 
graphic  coefficients  ordered  by  magnetic,  rather  than  geographic,  latitude.  The 
success  of  this  experimental,  data  based,  or  empirical  model,  in  representing  the 
actual  world-wide  foF2  is  due  to  the  large  amount  of  data  available  from  ionosondes 
in  many  regions  of  the  world.  Other  characteristics  of  this  model  are  discussed 
in  Dandekar.  *  * 

For  the  TEC  parameter,  data  availability  have  been,  and  will  likely  continue 
to  be,  much  more  sparse.  First,  TEC  measurements  have  been  made  generally  by 
using  measurements  of  Faraday  polarization  rotation  using  VHF  signals  of  oppor¬ 
tunity  transmitted  from  geostationary  satellite  telemetry  transmitters.  A  few  lunar 
reflected  Faraday  rotation  measurements  in  the  late  1950's  and  early  1960's  and 
the  TEC  obtained  from  a  few  low  orbit  satellites  did  not  contribute  significantly  to 
our  knowledge  of  world-wide  TEC  behavior,  at  least  not  for  modelling  average 
ionospheric  conditions.  Only  since  the  early  to  mid- 1960's  have  TEC  values  been 
obtained  on  a  more-or-less  regular  basis.  Even  today  fewer  than  one  dozen  stations 
regularly  contribute  TEC  data,  which  can  be  used  in  TEC  modelling  purposes,  to  a 
world  data  center. 

Fortunately,  most  of  the  contribution  to  TEC  comes  from  near  the  F2  region 
density  peak  where  models  of  foF2  are  available.  These  foF2  models  can  be  com¬ 
bined  with  some  limited  knowledge  of  topside  ionospheric  thickness  obtained  from 
topside  sounders,  along  with  topside  in-situ  density  measurements,  to  produce  a 

complete  ionospheric  height  profile  model.  The  most  well  known  of  these  models 
12 

is  the  one  by  Bent.  which  uses  ITS-78  coefficients  for  foF2,  and  topside  expo¬ 
nentials  for  computing  TEC.  A  representation  of  world-wide  average  behavior  of 
TEC  is  illustrated  in  Figure  6  for  2000  hrs  UT.  To  first  order  the  TEC  contours 
shown  in  Figure  6  move  westward  along  magnetic,  rather  than  geographic,  latitude 
lines,  at  the  earth’s  rotation  rate.  The  Bent  model  was  constructed  using  solar 
maximum  data  from  the  1968-69  period  and  had  to  be  adjusted  upward  somewhat  to 
account  for  the  much  higher  1979-1980  solar  maximum  than  that  of  1968-1969. 

This  adjustment  was  necessary  to  adequately  represent  the  actual  TEC  values  from 
stations  making  observations  in  March  1980,  which  was  near  the  maximum  of  the 
second  highest  solar  cycle  ever  recorded  in  the  over  200-year  history  of  solar 

10.  Barghausen,  A.  L. ,  Finner,  J.  W. ,  Proctor,  L.  L.  ,  and  Schultz,  L.  D.  (1969) 

Predicting  Longterm  Operational  Parameters  of  High-frequency  Skv-wave 

Telecommunication  Systems,  ESSA  Tech  Report,  ERL  110-ITS  78. 

11.  Dandekar,  B.  S.  ( 1982)  Ionospheric  Modeling.  AFGL-TR-82-0024,  ADA115243. 

12.  Llewellyn,  S.K. ,  and  Bent,  R.  B.  (1973)  Documentation  and  Description  of  the 

Bent  Ionospheric  Model,  AFC  RL-TR-73-0S57,  AD  772733. 
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cycle  observations.  The  Bent  model,  appropriately  adjusted  for  high  solar  cycle 

values,  does  however,  represent  fairly  well  the  average  behavior  of  TEC  for  many 

locations  tested.  Other  world-wide  ionospheric  electron  density  profile  models  from 

which  average  TEC  can  be  obtained  includes  one  by  Ching  and  Chiu,  ^  and  Chiu,  ^ 

Kohnlein,  ^  the  4-D  model,  ^  and  the  International  Reference  Ionosphere  (IRI).  ^ 

The  characteristics  of  some  of  these  models  are  described  in  Dandekar.  * 

Other  empirical  models  of  TEC  have  been  developed  directly  from  TEC  data 

alone,  though  these  have  necessarily  been  limited  in  temporal  and  geographic  extent 

by  the  available  data  base.  These  include  models  of  TEC  over  Europe  and  the 
1 8 

Mediterranean  for  low  and  medium  activity  portions  of  the  11-year  solar  cycle, 

and  a  model  of  TEC  over  the  Indian  sub-continent  for  both  solar  minimum  and  for 

1 9 

an  average  solar  maximum.  Models  of  the  slab  thickness  parameter,  the  ratio 

of  TEC / N  have  been  developed  for  specific  regions  such  as  the  one  for  northern 
max  90  21 

Europe  bv  Kersley,  and  one  for  the  eastern  USA  by  Klobuchar  and  Allen,  from 

which  TEC  can  be  obtained  from  a  model  of  foF2.  An  algorithm  designed  for  an 

approximate  50  percent  correction  to  world-wide  TEC,  for  use  in  an  advanced 

22 

navigation  system,  has  been  developed  by  Klobuchar. 

All  of  the  models  listed  here,  and  the  list  is  by  no  means  complete,  are  em¬ 
pirical  models  which  attempt  to  correct  for  average  TEC  behavior  only.  However, 
the  variability  from  average  TEC  behavior  can  be  large  and  may  be  important  to 
some  radio  wave  systems  which  must  propagate  signals  through  the  ionosphere. 

2.2  Temporal  Variability  of  TEC 

2.2.  1  VARIABILITY  FROM  MONTHLY  MEAN  TEC  VALUES 

The  ionosphere  is  a  weakly  ionized  plasma  and  the  resultant  TEC  is  a  function 
of  many  variables  including  solar  ionizing  radiation,  neutral  wind  and  electric  field 
effects,  neutral  composition  and  temperature  changes.  A  monthly  overplot  of  curves 
of  diurnal  changes  in  TEC  for  a  northern  mid-latitude  station  for  twelve  months 
during  a  solar  maximum  period  is  shown  in  Figure  7.  The  standard  deviation  from 
monthly  mean  diurnal  behavior  is  approximately  20-25  percent,  during  the  daytime 
hours  when  the  absolute  TEC  values  are  greatest.  Figure  8  shows  the  standard 
deviation  from  monthly  average  TEC  behavior  for  the  mid-day  hours  for  a  number 
of  stations  during  the  solar  maximum  period  1968-1969.  Again  20-25  percent  is 
a  good  value  for  the  standard  deviation  from  the  monthly  average  behavior.  The 
standard  deviation  is  somewhat  higher  during  the  nighttime  hours,  but  the  absolute 
TEX.'  values  are  much  lower  during  those  periods. 


(Due  to  the  large  number  of  references  cited  above,  they  will  not  be  listed  here. 
Refer  to  References,  page  29.  ) 
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PERCENTAGE  STANDARD  DEVIATIONS 
12  -  16  LT  TEC 
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Figure  8.  Percentage  Standard  Deviations  for  Daytime  TEC  From  the 
Stations  Indicated 


If  a  satellite  ranging  system  has  error  requirements  such  that  it  must  correct 

for  monthly  average  ionospheric  time  delay,  but  still  can  tolerate  the  approximate 

20-25  percent  variability  of  TEC  from  monthly  average  conditions,  approximately 

70-80  percent  of  the  ionospheric  effect  on  the  system  can  be  eliminated  by  the  use 

12 

of  an  average  TEC  model  such  as  the  one  constructed  by  Bent.  If  the  system  only 

requires  an  approximate  50  percent  rms  correction  of  the  ionospheric  time  delay 

22 

the  algorithm  developed  by  Klobuchar  can  be  used.  On  the  other  hand,  if  correc¬ 
tions  for  some  portion  of  the  remainder  of  the  ionospheric  time  delay  are  required, 
after  a  state  of  the  art  TEC  model,  such  as  the  Bent  one,  has  been  used  to  take 
out  the  monthly  mean  TEC,  then  the  short  term  (a  few  hours)  temporal  variability 
as  well  as  the  geographic  variability,  of  TEC  must  be  considered. 

2.2.2  SHORT  TERM  TEMPORAL  VARIABILITY  OF  TEC 

The  correlation  time  of  departures  of  TEC  from  monthly  average  curves  has 

23 

been  studied  by  Donateili  and  Allen.  They  concluded,  for  the  mid-latitude  station 
they  studied,  that  the  useful  prediction  time  was  a  function  of  local  time,  season, 
and  long  term  sunspot  activity.  However,  in  most  cases  they  found  no  significant 
improvement  over  the  use  of  monthly  mean  predictions  when  they  used  actual  data 
over  three  hours  old.  The  longest  useful  prediction  interval  occurred,  fortunately. 


23.  Donateili,  D,  C. ,  and  Allen,  R.  S.  ( 198 1 )  Time  cells  for  adaptive  predictions  of 

total  electron  content.  Radio  Sci.  16(No.  2) : 26 1  - 269. 

■■ 
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during  solar  maximum  daytime  hours  when  absolute  TEC  values  are  highest.  During 
solar  minimum  periods  their  useful  prediction  time  interval  was  often  as  short  as 
one  hour. 

23 

The  TEC  prediction  study  done  by  Donatelli  and  Allen  was  for  data  in  which 
the  predictions  were  for  the  same  geographic  location  and  direction  in  the  sky  as 
the  measurements.  If  the  prediction  is  for  a  different  location  the  temporal 
correlation  will  be  lower. 

2.2.3  GEOGRAPHIC  VARIABILITY  OF  TEC 

The  variability  of  TEC  at  the  same  local  time,  but  as  a  function  of  distance  has 

24 

been  studied  by  Klobuchar  and  Johanson.  They  used  TEC  data  from  two  sets  of 
stations,  one  aligned  approximately  along  an  east- west  direction,  with  the  other  set 
of  stations  aligned  along  an  approximate  north-south  direction.  Their  results  are 
shown  in  Figures  9a  and  9b  for  the  east-west  and  the  north-south  station  alignments, 
respectively.  No  significant  difference  in  correlation  distance  was  found  with 
season. 

The  percent  improvement,  P.  I. ,  in  TEC  from  the  average  value,  if  a  measured 
value  is  available  from  a  location  where  the  correlation  coefficient  between  the  two 
locations  is  r.  is  related  to  the  correlation  coefficient,  r,  by: 

P.  I.  =  100  *  (1  -  (1  -  r2)0' 5]  (see  Gautier  and  Zacharisen25)  . 

Note  that  a  correlation  coefficient  of  0.7  explains  only  29  percent  of  the  variance 
between  the  data  at  station  pairs;  hence  a  measurement  at  one  station,  with  a 
correlation  coefficient  of  0.7  between  data  sets  with  a  second  station  would  result 
in  an  improvement  at  the  second  station  over  the  average  predicted  value  of  only 
29  percent. 

2.3  TEC  in  the  Near-Equatorial  Region 

All  of  the  preceeding  sections  have  concentrated  on  the  behavior  of  TEC  in  the 
mid-latitude  regions  of  the  world,  mainly  because  most  of  the  available  data  are 
from  that  region.  The  near-equatorial  region  deserves  special  mention  due  to  the 
fact  that  the  highest  TEC  values  in  the  world  occur  in  this  region,  as  shown  in 
Figure  6.  This  region  extends  to  approximately  ±  20-25  degrees  on  either  of  the 
magnetic  equator,  with  the  highest  TEC  values  not  at  the  equator,  but  rather  at  the 

24.  Klobuchar,  J.A.,  and  Johanson,  J.  RI.  (1977)  Correlation  Distance  for  Mean 

Daytime  Electron  Content,  AFGL-TR-77-OlSS,  ADA048117. 

25.  Gautier,  T.  N. ,  and  Zacharisen,  D.  II.  (1965)  Use  of  space  and  time  correlations 

in  short  term  ionospheric  predictions,  IEEE  Conf,  Rec.  1st  Annual  IEEE 
Communications  Convention,  pp.  671-67?T 


23 


IONOSPHERIC  OISTANCE  (km) 

L  ATITUOE  SEPARATION  ■70*W  LONS  ITU  DC 


LONGITUDE  SEPARATION  AT  NORTHERN  MID  -  LATITUDES 


24 


Figure  9.  Correlation  Coefficient  vs  Station  Separation  in  Latitude  (9a)  and  Longitude  (9b) 


so  called  "equatorial  anomaly"  regions  located  at  approximately  ±  15  degrees  from 
the  magnetic  equator.  The  regions  of  highest  TEC  values  at  2000  hrs  UT  are 
rlearlv  seen  in  Figure  6  near  100  degrees  west  longitude. 

Most  of  the  day-to-day  geographic  variability  of  TEC  in  the  equatoriaL  anomaly 
region  during  solar  minimum  conditions  can  be  explained  by  the  variability  of 
equatorial  electrojet  strength.  Unfortunately,  no  such  similar  TEC  data  are 
available  for  solar  maximum. 

An  example  of  the  high  temporal  variability  of  TEC  for  solar  maximum  condi¬ 
tions  for  Ascension  Island,  a  station  located  near  the  peak  of  the  southern  TEC 
equatorial  crest  region,  is  shown  in  Figure  10.  Note  the  extremely  large- day-to- 
day  TEC  variability  in  the  afternoon  and  evening  hours  in  some  months.  Any  satel¬ 
lite  ranging  system  requiring  ionospheric  TEC  corrections  in  the  near-equatorial 
region  should  not  use  the  mid-latitude  standard  deviation  values  of  approximately 
20-25  percent  to  represent  the  variability  of  the  near-equatorial  region. 

2.4  TEC  in  the  Auroral  and  Polar  Cap  Regions 

Since  most  available  TEC  values  have  been  measured  from  radio  signals  trans¬ 
mitted  from  geostationary  satellites,  which  can  be  viewed  only  at  low  elevation 
angles  from  high  latitudes,  knowledge  of  the  variability  of  TEC  in  the  auroral  and 
polar  cap  regions  is  sparse.  In  the  American  longitude  sector,  where  the  magnetic 
latitudes  are  lowest  for  a  given  geographic  latitude,  there  is  considerable  TEC  data 
from  Goose  Bay,  Labrador.  At  Goose  Bay  the  aurora  passes  southward  even  during 
moderately  magnetically  disturbed  periods.  The  behavior  of  TEC  during  those 
periods  can  be  highly  irregular,  especially  during  the  nighttime  hours  when  TEC 
values  often  exhibit  rapid  changes  and  occasionally  even  exceed  the  daytime  maxi¬ 
mum  values  briefly.  While  the  occurrence  of  general  auroral  activity  may  be  pre¬ 
dictable,  the  specific  large  increases  in  TEC,  likely  due  to  auroral  precipitation, 
are  not  individually  predictable,  but  may  be  statistically  characterized  as  a  function 
of  magnetic  activity. 

In  the  polar  cap  region  a  negligible  amount  of  TEC  data  exists.  The  absolute 
TEC  values  are  likely  lower  in  this  region  than  in  the  mid-latitudes,  and  the 
variability  of  the  polar  cap  TEC  is  probably  very  high. 
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2.5  Protonospheric  Electron  Content 


Most  of  the  available  TEC  data  has  been  taken  by  measurements  of  Faraday 
rotation  of  single  frequency  radio  waves  transmitted  from  geostationary  satellites 
to  ground  observers.  The  electron  content  obtained  from  Faraday  rotation  observa¬ 
tions,  while  made  from  radio  waves  transmitted  from  satellites  at  geostationary 
satellite  height  above  the  earth's  surface,  includes  only  the  contribution  of  electrons 
up  to  heights  of  approximately  2000  km.  This  is  because  the  integrated  product  of 
the  longitudinal  component  of  the  earth's  magnetic  field  times  the  electron  density, 
above  approximately  2000  km,  is  negligible.  The  only  measurements  of  the  addi¬ 
tional  contribution  of  electrons  above  the  Faraday  maximum  height  have  been  made 

using  signals  from  an  ionospheric  beacon  on  the  geostationary  satellite  ATS-6. 

2  6 

Davies  has  reviewed  the  overall  results  of  the  ATS-6  experiment.  A  summary  of 

typical  protonospheric  electron  content  data  is  shown  in  Figure  11  taken  from 
27 

Klobuchar  et  al.  Note  that  the  protonospheric  values  are  fairly  low  in  absolute 
value. 

During  the  nighttime  hours  when  the  ionospheric  TEC  is  low,  the  protonospheric 
contribution  may  become  a  fairly  large  percentage  of  the  total  number  of  electrons 
between  a  satellite  at  geostationary  height  and  an  observer  on,  or  near  the  earth's 
surface.  Unfortunately,  no  protonospheric  electron  content  are  available  during 
solar  maximum  conditions. 

2.6  Short  Term  Variations  in  TEC 

The  time  rate  of  change  of  TEC,  in  addition  to  the  normal  diurnal  variations, 
also  has  periodic  variations  due  to  perturbations  of  the  ionospheric  F  region  from 
various  potential  Sources  as  geomagnetic  substorms,  meteorological  sources  such 
as  weather  fronts,  shock  waves  from  supersonic  aircraft,  volcanic  explosions, 
rocket  launches,  and  other  miscellaneous  sources.  While  these  short  term  varia¬ 
tions  in  TEC  cover  a  large  range  of  periods  and  amplitudes,  common  periods  range 
from  20  to  over  100  min  with  amplitudes  of  a  few  percent  of  the  background  TEC. 

A  10  percent  ionospheric  disturbance  with  respect  to  the  background  TEC  is  very 
*  28  29 

rare,  while  a  l  percent  TEC  perturbation  is  common.  Titheridge  and  Yeh  have 

made  studies  of  the  statistics  of  travelling  ionospheric  disturbances  (TID's),  in  TEC 

for  mid-latitude  regions. 

A  system  which  requires  correction  for  the  rate  of  change  of  TEC  cannot  rely 
on  models  of  TEC  to  provide  reliable  information  on  short  term  rate  of  change  of 
TEC  information,  and  can  use  available  TID  information  only  in  a  statistical  manner. 


(Due  to  the  large  number  of  references  cited  above,  they  will  not  be  listed  here. 
See  References,  page  29.) 
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The  only  recourse  for  a  system  significantly  affected  by  rate  of  change  of  TEC  is 
to  use  a  dual  frequency  measurement  technique  to  measure  directly  the  ionospheric 
contribution  to  range  rate. 
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Figure  11.  Monthly  Average  Plasmaspheric  Electron  Content  vs  Local 
Time  for  Aberystwyth,  Wales  (dashed  line)  and  for  Hamilton,  Massa¬ 
chusetts  (solid  line).  Also  ploted  are  values  for  Kiruna,  Sweden  for 
October  1975  and  from  Sao  Paulo,  Brazil  for  May  1975 


3.  CONCLUSIONS 

There  are  at  least  three  categories  of  systems  potentially  affected  by  iono¬ 
spheric  time  delay.  For  the  first  category  of  user  the  potential  systems  effects  may 
*  be  small,  at  least  under  any  naturally  occurring  worst  case  ionospheric  conditions. 

In  the  second  category,  a  user  may  require  a  nominal  correction  for  average 
ionospheric  time  delay,  but  is  able  to  tolerate  the  20-25  percent  standard  deviation 
.  from  average  conditions.  He  should  expect  at  least  a  50  percent  correction  for 
ionospheric  time  delay  effects  using  a  relatively  simple  time  delay  algorithm,  and 
up  to  70-80  percent  for  a  state  of  the  art,  fairly  complex  model,'  These  model 
corrections  can  be  improved  by  the  use  of  actual  ionospheric  measurements  within 
a  reasonable  temporal  and  spatial  frame.  For  the  third  category  of  user,  ionospheric 
model  corrections,  even  updated  with  near- real-time  measurements,  may  not  be 


28 


sufficient  to  correct  for  ionospheric  time  delay,  and  the  system  must  then  make  its 
own  ionospheric  correction.  Fortunately,  the  ionosphere  is  a  dispersive  medium 
and  the  use  of  identical  modulation  on  two,  widely-spaced  frequencies  will  allow  a 
direct  measurement  to  be  made  of  ionospheric  range  delay.  Two  coherently-derived 
carrier  frequencies  may  be  used  to  obtain  accurate  time  rate  of  change  information 
for  TEC.  Details  of  measuring  ionospheric  effects  directly  by  a  system's  use  of 
multiple  frequencies  are  available  in  Burns  and  Fremouw.  4 
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